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Abstract 

We present exact expressions for the Sagnac effect of Godel's Universe. For this 
purpose we first derive a formula for the Sagnac time delay along a circular path in 
the presence of an arbitrary stationary metric in cylindrical coordinates. We then 
apply this result to Godel's metric for two different experimental situations: First, 
the light source and the detector are at rest relative to the matter generating the 
gravitational field. In this case we find an expression that is formally equivalent 
to the familiar nonrelativistic Sagnac time delay. Second, the light source and 
the detector are rotating relative to the matter. Here we show that for a special 
rotation rate of the detector the Sagnac time delay vanishes. Finally we propose a 
formulation of the Sagnac time delay in terms of invariant physical quantities. We 
show that this result is very close to the analogous formula of the Sagnac time delay 
of a rotating coordinate system in Minkowski spacetime. 



* email : endre . kaj ari@physik. uni-ulm. de 



1 



1 INTRODUCTION 



2 



1 Introduction 

Phenomena in rotating coordinate systems have fascinated scientists since hundreds of 
years. Three examples may serve as an illustration of this statement: i) The Coriolis 
force is instrumental in the demonstration of the rotation of the earth with the help of 
the Foucault pendulum pj. ii) The Sagnac effect [21 El IH El EI measured by the Michelson- 
Gale-interfero meter [7j is the analogous optical tool, in) Mach's principle 0131101 ushers 
in a fresh view on the relativity of rotation. In the present paper we combine these three 
concepts and calculate the Sagnac effect of Godel's Universe. 

Einstein's theory of General Relativity JT] predicts gravito-magnetic forces ^21 ^1 ITij 
due to the rotation of massive objects. These forces give rise to the precession of nodal 
lines ^Hl of two orbiting LAGEOS-satellites proving the first measurement of the Lense- 
Thirring effect [THl IT71 ITHl IT^ . The newly emerging field of atom optics has opened a new 
arena (201 for measuring the dipole character of the Lense-Thirring field with the help of 
an atom gyroscope |2I1 1221 123 121 • 

Rotation is not limited to a coordinate system or individual masses but even the universe 
can display features of rotation. Indeed, in 1949 Kurt Godel |2S1 123 12Z| derived an exact 
solution of Einstein's field equations, in which a homogeneous mass distribution rotates 
around every point in space. This solution shows rather unusual properties |2H] such as 
closed time like world-lines. 

In a recent paper [23 we have evaluated the Sagnac effect of Godel's Universe measured 
in a laboratory size interferometer. For this purpose we have used the linearized Godel 
metric, which can be approximated by the flat spacetime metric of a rotating coordinate 
system. In the present paper we calculate the exact Sagnac effect for the Godel metric. 

1.1 Why Godel's Metric? 

Godel's Universe has many fascinating features. Indeed, the inherent rotation of this 
Universe is only one interesting aspect. Even more intriguing is the lack of a global time 
ordering and the existence of closed time like world lines giving rise to the possibility of 
time travel. Such causal problems emerge also in other exact solutions of Einstein's field 
equations, such as the Kerr metric. However, Godel's metric has the advantage that it is 
of rather compact form and most calculations can be carried out analytically. 

It is commonly accepted that this model does not provide a matching description of our 
observed universe. Nevertheless, there is the emerging field of experimental cosmology 
in the laboratory jSU]. In particular, the examination of wave phenomena in curved 
spacetimes is a focus of research. For example, optical analogues of black holes have been 
proposed by studying light propagation in moving media [SIl or sound propagation in 
condensed matter systems [221 E3]- In this context, it is important to ask how far one can 
extent such analogues in general. Godel's Universe can shed some light on this problem, 
since the existence of closed time like world lines curtails the expectation of a globally 
valid experimental analogue. Thus, it would be interesting to see, where and how the 
experimental realization ceases to exist. 
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Motivated by this idea we have initiated a study of wave propagation for the source free 
Maxwell equations in Godel's metric jS^ IHS]. However, in the present paper we will focus 
on the Sagnac time delay within the limit of geometrical optics and draw the comparison 
to the Sagnac time delay in a rotating frame in flat Minkowski spacetime jSHI- In order 
to keep the paper self-contained and in view of the fact that this issue brings together 
researches from atomic, molecular and optical physics with experts in General Relativity 
we have used this opportunity to combine our study of the Sagnac effect with a mini 
review of Godel's Universe. 

1.2 Outline of the Paper 

Our paper is organized as follows: In Section |21 we first briefly review the essential features 
of the Sagnac effect and then derive an exact expression for the Sagnac time delay for a 
time independent metric in cylindrical coordinates. We dedicate Section El to a discussion 
of the Godel metric. We then apply in Section 0] the expression of the Sagnac time delay 
to Godel's metric. In SectionElwe reformulate the Sagnac effect purely in terms of proper 
time delays of light pulses and compare it to the analogous expression for a rotating 
frame in flat Minkowski spacetime. We conclude by an outlook to possible laboratory 
realizations of the local light propagation in Godel's Universe. 

2 Sagnac Effect: Basics 

In the present section we first recall the familiar Sagnac effect and then derive an exact 
expression for the Sagnac time delay in terms of metric coefficients. We conclude by 
making contact with the familiar Sagnac effect using the metric of a rotating coordinate 
system in fiat Minkowski spacetime. 

2.1 Familiar Sagnac Time Delay 

In 1913 George Sagnac performed an experiment [21 IS], which he interpreted as an veri- 
fication of the existence of the ether. In order to bring out the essential features of this 
phenomenon we consider two counter-propagating light pulses which travel on a circle of 
radius due to an appropriate array of mirrors or a glass fiber. When the setup is at 
rest the two counter-propagating pulses arrive at the same time at the point of emission. 
However, when the arrangement rotates there is a time delay between them. This Sagnac 
time delay Ats follows from an elementary, classical argument and reads 

Ats ~ -^^rl (1) 

Hence, Ats is proportional to the area nr^ enclosed by the light beams and the rotation 
rate Qr of the ring with respect to the fiat Minkowski spacetime. Moreover, we note that 
this result is an approximation and higher order corrections in Qf^ will arise as discussed 
in Section 1231 
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The Sagnac effect is the basis of modem navigational systems. Moreover, it can also 
be generalized in the framework of General Relativity IHZI , as presented in the next 
section. In particular, the Sagnac effect is closely related to to the Hannay angle jHEI and 
the synchronization of clocks [SHI along a closed path. 



2.2 Time Delay for a Stationary Metric 

We now derive an exact expression for the Sagnac time delay Ats in terms of the time inde- 
pendent metric coefficients (^^j, expressed in cylindrical coordinates 

{t,r,(j),z). For a general treatment of the Sagnac effect of two counterpropagating light 
rays along any spatially closed path in an arbitrary metric we refer to ^40j. However, 
inhere we choose a particularly simple configuration and consider the propagation along 
a circle of radius tq in the z = zq plane resulting in dx^ = dr = and dx^ = dz = 0. For 
this situation the line element reads 

ds"^ = g^iydx'^dx" = goo{dx^Y + 2 dx^dx^ + g22idx'^Y, (2) 

and vanishes for light ds"^ = 0. 
The angular velocity 



dx'^ , , 

of the light beam follows from the quadratic equation 

gnLO^ + 25(020; + 5'oo = 0, 

giving rise to the two velocities 

= ± Jf^V-^. (4) 

5'22 V Vfi'22/ g22 

In order to get two solutions, which correspond to two ordinary counter-propagating light 
beams, we restrict ourselves to spacetime regions G, in which the above metric coefficients 
satisfy the conditions 

^700 > 0, (722 < V G G. (5) 

These conditions immediately imply that the angular velocity a;+ is positive and cu. 
is negative for all events in G, so that the monotony of the corresponding solutions 
x\[x'^^ is guaranteed. The condition ^fQO > is furthermore important to allow for an 
observer resting relative to the chosen spatial coordinates and measuring the proper time, 
as discussed e.g. in |i4T| . 

We can find the coordinate times x\ of the first return of the counter-propagating light 
rays to their starting point x^ = 0o by integrating equation Q over one period 2tt in 
positive and negative angular direction, respectively. Hence, we arrive at 

<j>0±2TT 2tt 

xl= I f^ = ±/^, (6) 

J UJ± J UJ± 

<bn 
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where in the second step we have made use of the periodicity of the metric coefficients in 
the coordinate x^. 

When we recall the connection dr = a/5'oo('^o, 0o, ^o) (ia;°/c between the coordinate time x° 
and the proper time r measured by an observer resting at (ro, ^o, zq), the corresponding 
proper times t± of the incoming light rays read 



2n 



C J U}± 



(7) 



The positive sign in r+ corresponds to the light pulse propagating in the positive angular 
direction, whereas the negative sign in r_ denotes the proper time of the light pulse 
traveling in the negative angular direction. 

The Sagnac proper time delay Ats follows as 

Ats = (r+ - r_) = -\^goo{ro,4>o,Zo) / — — — dx^, 



which reduces with help of the exphcit expressions dH) for uj± to 

2tt 



Ats = --Vgoo{ro,4)o,zo) [ — dx^. 
c J goo 



If we further assume, that the metric coefficients (702 and g^o do not depend on the angular 
coordinate x^, then the expression (jHl) for the Sagnac time delay reduces to the compact 
formula 

Ar. = -1^ ^. (9) 



For negative values of Ats the light pulse which propagates in the positive angular direc- 
tion returns before the other pulse and vice versa for positive Ats- 

We conclude by noting that the Sagnac time delay is by construction only form- 
invariant under the special class of coordinate transformations 

'0 _ '0( 1 3\ n _ '1(^1 ™3n /2 _ 2 /3 _ '3/1 3\ 

which change neither the frame of reference nor the angular coordinate. 



2.3 Time Delay in a Rotating Frame 

We now want to apply formula (jH)) of the Sagnac time delay to the metric of a rotat- 
ing coordinate frame in flat Minkowski spacetime. The line element in flat Minkowski 
spacetime reads in cylindrical coordinates = {t, r, 0, z) 



ds^ = c^dt^ - dr"^ - r'^dcj? - dz^. 



(10) 
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The coordinates x''^ = {t', r', 0', z') of a reference frame rotating with a rate Qr > are 
defined by the transformation equations 

t = t', r = r', (p = (p' + njit', z = z', 

and give rise to a hne element ds"^ = g'^^dx'^dx"^ of the form 

ds^ = (c^ - r'^n\) dt'^ - dr'^ - r'^d<f)'^ - dz'^ - 2r'^QR dt'dcj)'. (11) 

To satisfy the conditions © we have to restrict ourselves to the spacetime region 

G" = jo < r' < < 0' < 27r, -oo < t\ z < ooj . (12) 

When we now assume that the counter-propagating light rays travel along a circle with 
radius Tq and substitute the metric coefficients q'qq and §'^2 into Q, the Sagnac time delay 
in the rotating frame reads 



In the first approximation of {r'Vtii/c) this expression reduces to the time delay of the 
familiar Sagnac effect dH). 

However, also the limit of large rotation rates, where the square root in the denomi- 
nator gets important is of interest. Indeed, in Section El we present an expression for 
the Sagnac time delay in Godel's Universe, which on first sight looks very different from 
(|13|) . Nevertheless, a closer view reveals that this formula is numerically very close to the 
Minkowskian Sagnac time delay. 



3 Essential Features of Godel's Universe 

In order to gain some insight into the intricacies of Godel's Universe we briefly review 
two representations of Godel's metric that are convenient for our analysis of the Sagnac 
effect. Moreover, we sketch the conditions under which this metric solves Einstein's fleld 
equations. We then present tensorial quantities characterizing the time like velocity fleld 
of the matter generating the Godel metric. Furthermore, we make contact with the 
metric of flat spacetime in cylindrical coordinates and mention the symmetries and the 
causal structure of Godel's metric. We conclude this section by summarizing special null 
geodesies, which will be important in expressing the Sagnac time delay by measurable 
quantities. 

3.1 Line Element and Einstein's Field Equations 

The line element ds^ = g^^, dx^dx'^ given by Godel |25j in 1949 in dimensionless, cylindrical 
coordinates x'^ = (x°, x^, x^, x^) = (t, f, 0, z) with 

Gg = {-oo < t < oo, < f < oo, < < 27r, -oo <z< oo} , (14) 
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dP — dr"^ — (sinh^ r — sinh*^ r) dcj)^ — dz^ + 2v^ sinh^ ; 



^dt. 



(15) 



has the form 

ds^ 

The parameter a > has the unit of a length. 

In the next section we derive a simple expression for the Sagnac effect. For this purpose, 
it is convenient to use a slightly different form of Godel's metric, which can be obtained 
from (fT3j) by the coordinate transformation 



2a- ^ . , - 

t = — t, r = 2a smhr, 



z = 2az. 



(16) 



The resulting line element reads 

dr'^ 



ds' = c^dt^ 



^ \ 2 

2a 



- dz^ + 2r^ 



These new coordinates have now physical dimensions. 
We can convince ourselves that the metric coefficients 



{9, 



\ 



C 



.2 c 









„2 c 



V2a 





-r^ ( 1 



f-V 

\2a/ 



V2c 



-dt ( 







corresponding to the line element (fTTj) indeed solve Einstein's field equations 



0\ 





-1 



(17) 



by calculating explicitly the Ricci tensor R^^, and the scalar curvature i?, defined in Ap- 
pendix ^ Here A denotes the cosmological constant and /t = (87rG)/(c^) with Newton's 
gravitational constant G and the velocity of light c. 

The energy momentum tensor T^i, of an ideal fluid with mass density p and pressure p 
reads 

T^,u = (p + ^) u^u^ - pg^^ . 
With the four-velocity = (1, 0, 0, 0) of the matter generating the field, we find 



R 



9 9 * 



R 



When we substitute these expressions into the field equations, we arrive at the two rela- 
tions 

K,p = A + — — , 



9 9 ' 



2a? 



which couple the length scale a to the density p, the pressure p and the cosmological 
constant A. 
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3.2 Time Like Velocity Field of the Ideal Fluid 

The essential properties of the motion of the ideal fluid generating the field are char- 
acterized by the tensorial quantities 9, and Uafi representing the volume expansion, 
the shear tensor and the rotation tensor, respectively [2H1 From the four- velocity 
= (1, 0, 0, 0) and the acceleration = u^.^u'^ = of the congruence of time like curves 
belonging to the ideal fluid we find a vanishing volume expansion 



and a vanishing shear tensor 



Here we have introduced the projection tensor 

together with = {u^-y + u^-^) /2, where the semicolon denotes the covariant deriva- 
tive. 

However, we arrive at a non-vanishing rotation tensor 

u^p ^ P^^P^p uy,,,^ = = r-|- [5l5l - 5l5l) . (19) 

Here the square brackets are defined by U[^-y] = {u^-^u — Uu-^^) /2 and the comma denotes 

the partial derivative. 

The corresponding rotation vector 

u;'^^ls-^^'u,u,.,s = ^6- (20) 

and the rotation scalar 

no ^ = ^ > °- ^^^^ 

are constant in every point of Godel's Universe. 

In the limit a — oo all rotation quantities (jl9|) . (j2(jp and (|2ip vanish. Moreover, for 
r/(2a) <^ 1 we find that the line element (fTTjl can be approximated in first order by 

ds^ = c^dt^ - dr^ - r^rf^^ - dz^ + 2r'^VLGdt d(f) + 0(^1). 

While the zeroth order approximation corresponds to the line element fllUp of flat 
spacetime in cylindrical coordinates the linear correction is reminiscent of the line element 
fjll|) of a rotating coordinate frame in Minkowski spacetime. 
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3.3 Killing Vectors and Symmetry 

If a spacetime manifold with the metric Qf^i, possesses symmetries, then they can be 
characterized by a special class of coordinate transformations x'" = x'°'{x^) which satisfy 
the condition 

9',.ix'-) = 9,.ix'n- (22) 

Such a coordinate transformation is called isometry. In particular, the infinitesimal isome- 
tries 

x"^ = x'^ + eCix^) (23) 

with e <^ 1 are of special interest, since every continuous isometry can be constructed 
successively by these infinitesimal isometries. 

The infinitesimal transformation (j23j) together with (j22p and the transformation law of a 
metric yields the condition 

Up + ^(3;a = (24) 

for the Killing vector field ^"(x^). The solutions of this linear system of partial differ- 
ential equations characterize the symmetry of a given metric. 

In the case of Godel's metric we find five Killing vectors as solutions of the Killing equa- 
tion ()24|) . Three of them are immediately found from ()22j). since Godel's metric does 
not depend explicitly on the coordinates (t,(f),z). With the constants A, B,C, D, E the 
complete solution of reads 

C{t, r, 0, z) = A5^ + B5^ + C5^ + Ddr, 0) + EC (r, - f ) , (25) 

with 

+ (^)') 

/ 

Since ()25|) contains a time like Killing vector Godel's Universe is stationary. Moreover it 
is also spatially homogeneous. However, a more detailed analysis of the Killing vectors 
shows, that Godel's metric is not static and also not isotropic. The latter feature is due 
to the existence of a rotational axis giving rise to a rotational symmetry in the z = const. 
planes. 

3.4 Causal Structure 

The non-vanishing rotation scalar has a dramatic consequence for the causal structure of 
Godel's Universe. In order to gain some insight into this feature, it is useful to consider 
infinitesimal light cones at different spatial points. Figure Q depicts such an arrangement. 
The cylindrical coordinates {t, r, 0, z) are embedded for illustration in a Cartesian frame 



e 



(t,x = r cos (j),y = r sin 0, z) 



(26) 
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Figure 1: Light cones in Godel's metric represented in the z = plane. The middle circle 
of critical radius = 2a separates the domains of different causal behavior. At every 
point of the inner domain the light cones lie outside of the t = plane as exemplified by 
the cones along the inner circle. In contrast, outside of the critical circle the cones cut 
through this plane. 



and the third spatial coordinate z is suppressed in the figure. For every instant of time 
there exists a domain in space with closed time like world-lines. To visualize this property, 
we consider the orientation of some selected light cones with respect to the plane of 
constant coordinate time t. 

At the origin of the coordinate system the axis of the light cone is orthogonal to the plane. 
As we move away from the origin the light cones start to tilt, as indicated in Figure ^ 
At the critical Godel radius = 2a, represented by the middle circle, the light cones are 
tangential to the plane of constant coordinate time t. This circle of radius is a light like 
curve. Outside this critical radius the inclination of the light cones increases further and 
allows the existence of closed time like curves, as shown by the outer circle in Figure ^ 
It is this peculiar feature of the causal structure which permits to connect two arbitrary 
events of spacetime by a time like curve, irrespectively of their ordering in the chosen 
coordinate time t. Indeed, we can start from a point within the inner circle and cross 
the critical Godel radius to explore the world beyond this border. During this trespassing 
we take a time like world line which spirals downwards into the past. Having regressed 
long enough on this trajectory we can finally return to the original causal domain, thus 
arriving before departing. 

We conclude by noting that closed time like world-lines are not limited to Godel's Universe. 
They also appear in other exact solutions of Einstein's field equations for rotating mass 
distributions. Examples include the rotating Kerr black hole j44lH3] and the van Stockum 
rotating dust cylinder Most recently is has been shown |17], that the gravitational 
field of a solenoid of light, that is a light beam bent along helical path also exhibits such 
closed time like world-lines. 
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3.5 Null Geodesies 

In the presence of a metric, the free motion of particles or the propagation of hght rays 
is described by the geodesic equations 

d'^x'^ „„ dx"" dx^ 



+ r'l«-Tr-jr = <X = o (27) 



and by the condition 



dX^ d\ dX 



For the definition of the Christoffel symbols F'^^ we refer to Appendix ^ and x^(A) = 
(t(A), r(A), 0(A), ^(A)) is the path of the particle or light beam. For massive particles 
e denotes the velocity c of light and the curve parameter A represents its proper time. 
Moreover, the tangent vector = dx^/dX is the four- velocity of the particle. In the case 
of a light beam we have to set e = and the curve parameter A has no physical meaning. 

The solutions of the geodesic equations for Godel's metric were first given by W. Kundt 
|48j . independently examined by S. Chandrasekhar and J. P. Wright [IHl and discussed 
in detail by Novello et. al. [50j. 

In Section El we express the Sagnac time delay in measurable quantities only. For this 
purpose we need to know the geodesic motion of light for special initial conditions. Indeed, 
the light pulse shall start for the initial curve parameter Aq = at the point r(0) = 
0, ^(0) = 0. Furthermore, the light ray shall have a vanishing 2;-component of the initial 
velocity, that is ■u^(O) = 0. Since we start at the origin, the radial velocity has to be 
positive, that is u^{0) > 0. 

In Appendix |B] we outline a procedure for obtaining the general solution of the geodesic 
equations ^T7\i for Godel's metric and integrate these equations subjected to these initial 
conditions. 

For the radial coordinate we find the expression 

"^^^-'-(|r^A)|, (28) 



2a 

where we have introduced the abbreviation 



sm 



V2u\0)nG. (29) 



Hence, when the curve parameter A reaches the value Ac = {2Ti)/ri the light pulse returns 
again to the origin r = 0. 

The coordinate time along these geodesies reads 

t(A) = -M°(0) A + [arctan (^2 tan {r]X/2)) + m(A) tt] , (30) 
where the integer 



m(A) 



r/A 1 
2^ 2 



(31) 



represents the greatest integer less than or equal to the number inside the brackets. 
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Within one cycle, < A < 271/1] , the coordinate time t{X) increases by the time interval 

V27r 



Atr 



G 



(^-1). 



For the angle coordinate 0(A) the integration yields 

0(A) = 0(0) + arctan (^^2 tan {r]X/2)j + [m(A) - m(A - vr/r^)] vr. (32) 

We conclude this section by illustrating the null geodesies for the special initial conditions 
in Figure El The left side displays the geodesic motion in the (x, ?/)-plane embedded in 



t/At, 





a b 



Figure 2: Null geodesies in the (x,?/)-plane (left) and coordinate time t{X) extending 
over 2 cycles (right) for the special initial conditions r(0) = 0, u^{0) = and 0(0) = 
(0, 27r/3, 47r/3). The three values a, b, c of the curve parameter A mark the three different 
positions on a null geodesic on the left. 



the coordinate frame ()2(ij) . The light signals emitted at the origin r = cycle in the 
positive angular direction. The right side of the figure shows the coordinate time t{X) for 
the increasing curve parameter A. The inflection points in the time coordinate appear 
when the light ray touches the critical Godel radius. 



4 Sagnac Time Delay in Godel's Universe 

The goal of the present section is to derive an exact expression for the Sagnac time delay 
of two counter-propagating light rays on a circle in the presence of Godel's metric. We 
also analyze the Sagnac effect observed by a detector rotating relative to the ideal fluid. 
Finally we briefly discuss the dependence of the Sagnac time delay on the choice of the 
spatial coordinates. 

4.1 Sagnac Effect in the Rest Frame of the Ideal Fluid 

We start by considering a situation in which the emitter and the detector are at rest in 
the coordinate frame of the ideal fluid. The light pulses propagate along a circle of radius 
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r = ro in the plane z = zq. Recalling Godel's metric (fTHj) and the rotation scalar (PT|) the 
relevant metric coefficients read 

5-00 = c^ 5-02 = rl^G. (33) 

In order to satisfy the conditions (0), we have to restrict ourselves to the spacetime region 

G = {0 < r < 2a, < < 27r, -oo < t, z < oo} . 

If we do not confine ourselves to this domain the coordinate time needed by one of the 
light pulses to return to the starting point will be negative for radii larger than the critical 
Godel radius re = 2a. This feature follows from the orientation of the future light cones 
in Figure H 

Furthermore an observer resting with the detector at tq > tq would measure a negative 
proper time of the returning light ray, which travels backward in the coordinate time. By 
restricting our experimental setup to the region G we avoid such alien situations in which 
the light pulse returns before it is emitted with respect to the chosen coordinate time. 
When we substitute the coefficients flS^ into formula the Sagnac time delay for Godel's 
metric reads 

Ar5 = -^T^rl (34) 

Since this expression is negative, we conclude from the definition of the proper time Arg, 
that the light pulse propagating in the positive angular direction will always return to its 
starting point before the pulse propagating in the negative direction. 

It is instructive to compare this result to the familiar Sagnac time delay At5 given by 
(ID). On ffist sight the absolute value of the Sagnac time delay (jH^ in Godel's Universe 
seems to be identical to the familiar Sagnac effect (0). However, when we recall that the 
coordinate tq in Godel's Universe does not represent a proper distance, we recognize that 
the similarity only arises from the special choice of our spatial coordinates. We will return 
to this point in Section lOl 

We conclude by noting that the gravitational time delay in Godel's Universe has also 
been calculated in [21]. However, the expression given in that paper is quadratic in the 
rotation scalar. 

4.2 Sagnac Effect Measured by a Rotating Detector 

We now analyze the Sagnac time delay for a slightly different experimental situation. 
The light source emitting the two counter-propagating light rays and the detector are 
now no longer at rest relative to the ideal fluid, but rotate relative to it. We denote the 
corresponding rotation rate by Due to this additional rotation the Sagnac time delay 
will also depend on fi^). After providing an explicit expression for the Sagnac time delay 
for this situation we choose fi^i such, that this time delay vanishes exactly. 

We start from the metric (jl8|) and perform the coordinate transformation 



t = t' , r = r', = 0' + Qot', z = z' 



4 SAGNAC TIME DELAY IN GODEL 'S UNIVERSE 



14 



into the rotating frame of the detector. In the new coordinates the hne element takes the 
form 



4aV " -) 1 + [Q 



1 - ^ J d(P'^ - dz'^ + 2r'2 -[I- ^ J ^dJ dt' d^' . 

The conditions (0) for the corresponding metric coefficients in the rotating frame lead to 
the restrictions 



and 

< r' < re 

on the radial coordinate r' and the rotation rate Vtr,. Indeed, these conditions lead to a 
nontrivial region G of allowed radii r' with respect to the chosen Vtr,. 

We again denote the radius of the circular light path by Tq. When we insert the relevant 
coefficients q'qq and 5^02 into formula © and take into account the rotation scalar (PT|) we 
find the Sagnac time delay 



Ar's 



The special choice 



^ -%T- (35) 



'0 

4a2 



of the rotation rate of the detector leads to a vanishing time delay Ar^ between the two 
counter-propagating light pulses. Hence, for infinitesimal radii Tq, the rotation velocity 
Qd of the detector is equal to the rotation scalar of Godel's Universe. 



4.3 A Different Choice of Spatial Coordinates 

We now return to the Sagnac time delay measured in the rest-frame of the ideal fiuid of 
Section I^!T1 According to equation (j34p the Sagnac time delay Ats depends in a linear 
way on the rotation scalar of Godel's Universe. 

We emphasize that this linearity is only due to the special choice of our space coordinates 
(fTB|) in the line element (fTTj) . Indeed, we can also express the Sagnac time delay Ats 
in the dimensionless space coordinates of the original Godel line element fllSj) . where 
the coordinate transformation (fTB|) does not change the frame of reference. In these 
coordinates (fT^ the Sagnac time delay reads 

Ats = -jt- sinh^ tq. (36) 
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Although in this representation the dependence of Ats on differs from that of equation 
fl34|) both formulae contain the same physics. This invariance is hidden behind the explicit 
dependence of Ats in (I34|l and (j36|) on the spatial coordinates, which have no immediate 
physical meaning. In the next section we are going to rectify this problem by introducing 
an invariant formulation via an operational definition of the radial coordinate. 

5 Invariant Formulation of the Sagnac Time Delay 

The different formulae ()34|) and (jHUj) for the Sagnac time delay bring out most clearly the 
question: How can we reformulate the radial coordinate of the mirrors in a measurable 
quantity? 

In the present section we outline a measurement strategy to address this question and then 
reformulate the Sagnac time delay solely in terms of proper times of light signals. The 
obtained relation allows for an interesting comparison to expression (jl3|) for the Sagnac 
time delay in a rotating frame of reference. 

5.1 Operational Definition of the Radial Coordinate 

All expressions for the Sagnac time delay derived so far contain a radial coordinate. 
What is a measurement strategy for obtaining the radius tq of the circular path of the 
light pulses? 

We answer this question in the spirit of [TT|IH^ IHH]. We send light signals along their null 
geodesies from the center of the circle to every mirror M in the z = plane, as illustrated 
on the left side of Figure El The mirrors are arranged in such a way as to reflect the light 
back to the point of emission. Moreover they are lined up in the z = plane as to ensure 
that the pulses emitted simultaneously at the center and reflected from the mirrors all 
return to the center at the same time. 

Hence, the radius tq can be expressed by the proper time Atm between the emission and 
detection of the reflected pulses measured by an observer resting at the origin r = 0. 
Once the circular setup of mirrors has been established by this operational procedure the 
mirrors have to be readjusted to guide the counter-propagating light rays in this circular 
Sagnac configuration with radius ro, as indicated in the right of Figure El 
We conclude by emphasizing that our experimental implementation is only the simplest 
conceivable model. It is a gedanken experiment, which even requires an infinite amount of 
mirrors. However, a practical realization has many caveats which need to be considered, 
such as the finite number of mirrors, the propagation of the light pulses on null geodesies 
in-between and the applicability of geometrical optics. In fact, the case of a finite number 
of mirrors has been carried out in principle by ^01 [33] , where the Sagnac time delay in a 
rotating frame in Minkowski spacetime serves as an example of the given method. 
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M 
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Figure 3: Illustration of the measurement procedure for the radius tq of a circle (left) 
and of the Sagnac time delay Ats (right). Light rays start at the origin and propagate 
on null geodesies to the mirrors (M), where they are reflected back to the origin (left). 
The proper time Atm between emission and arrival of the light pulses is measured by an 
observer resting at the origin. The figure on the right illustrates the typical experimental 
arrangement for the Sagnac effect with a detector (D) on the circle. The straight lines 
between two consecutive mirrors are an approximation to the corresponding null geodesies. 
However, the curves shown on the left are exact null geodesies (see Appendix IB]) . 



5.2 Invariant Formulation 

We are now in a position to establish a connection between the Sagnac time delay Ats 
and the rotation scalar iii terms of the proper time interval Atm, measured by an 
observer resting at the origin r = of the coordinate system. We start by calculating 
the coordinate time for the light pulses to propagate to the mirrors and back. For this 
purpose we recall ()28|) - (jHUj) . We then translate this time interval into proper time and 
express the radius tq by the Sagnac time delay making use of equation (jMj) . 

Since our mirrors are positioned at the radius tq < tq the curve parameter Aq corre- 
sponding to this radius can be found from equation (|28|) . Within a period < r^A < 27r 
two values of A correspond to tq and we choose the first solution Aq within the interval 
< r^Ao < vr. The coordinate time between the emission of a light pulse at r = and the 
reflection of it from r = tq follows from ()30p and ()31|) and reads 



Since Godel's metric is stationary, spatially homogeneous and possesses rotational sym- 
metry in the planes z = const., it takes twice the time t(Ao) to travel to and from the 
mirrors. For an observer resting at the origin r = the proper time Atm between emission 
and arrival of the light pulses is given hj ds = cdr = cdt, and consequently 




(37) 



Atm = 2t(Ao). 
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When we use equation (^Hj) for Aq in (|^ we arrive at the expression 

VIgAtm = -2V2 arcsin (^^) + 4 arctan [ ^^"^ ^ ] (3^ 



for the proper time Atm- Here we have recalled the definition (jSH))- 

In order to make the connection to the absolute value of the Sagnac time delay Arg, we 

recall from (jH^ the relation 

keeping in mind, that the negative value of Ats is only due to a faster propagation of the 
light pulse in the positive angular direction. 

This formula allows us to express the radius tq in in terms of the absolute value \Ats\ 
of the Sagnac time delay. In order to avoid the appearance of the square root of Ats it 
is convenient to introduce the dimensionless parameters 

Al^lna\Ars\^{^y (39) 



A^, = ^^gAtm. (40) 



with < Aq < 1 and 



When we replace the proper time Atm and the radius tq in (jHaj) by their corresponding 
dimensionless parameters (j39p and (j4(Jj) we arrive at the transcendental equation 

47rA^ = 2 arctan ( ^^^"^ ) - ^2 arcsin As. (41) 
\y/l-AlJ 

This relation is quite a remarkable result since it provides an invariant formulation of the 
Sagnac time delay in Godel's Universe. Indeed, for a given value of the Am, that is for 
the radius of the circular path measured in propagation time Atm of light, the solution 
of this equation is the scaled Sagnac time delay A^. Since we have scaled all proper 
times in terms of the rotation scalar the transcendental equation contains no parameters 
of Godel's Universe. The information about the metric reflects itself solely in the form of 
the equation. 



5.3 Comparison to the Rotating Frame in Flat Spacetime 

We now compare the transcendental equation (|4H) for the Sagnac time delay of Godel's 
metric with the corresponding formula for the time delay (fT^ found in a rotating frame 
of reference in Minkowski spacetime. 

The first step consists of replacing the radial coordinate in equation p3p by a measurable 
time. For this purpose we apply the same measurement strategy to the rotating frame as 
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discussed in 15.11 We then substitute the radial coordinate Tq by the proper time interval 
At'^j of a light ray, which returns after reflection from the circular mirror arrangement at 
to the origin. Since the observer at the origin is at rest relative to the inertial frame of 
Minkowski spacetime, we find immediately the proper time 



_ 2r'n 

M — 



At' - ° 



c 

For the comparison of the Sagnac time delay ()13|) with the expression fl41|) it is useful to 
introduce the dimensionless parameters 



and 



A'^ = l^^R^r'^ (42) 



AZ ^ ^nnAr',,, (43) 



in complete analogy to the formulae (jH^ and (jin|) . 

Using these parameters we can cast the Sagnac time delay ()13|) into the form 



2A« 



(4vrA',2,) 



1 - (4vrA'|,) 



which in contrast to ()41|) is an explicit formula for A'^ in terms of A\_f. 

In order to obtain a formula analogous to (PT|) we solve this equation for AnA'^^, which 

yields 



47rA'^, = v^A'^y y^l + A'^ - A'|. (44) 

Due to our restriction (|T^ to the spacetime region G' the parameter A'^ has to satisfy 
the condition 

< 47rA'^_^ < 1. 

Figure m illustrates the differences in the curves defined by ()4ip and ()44|) . We find that 
for our measurement strategy the difference between the solid and dotted curve is very 
small. This is surprising, as the underlying physical systems with Godel's Universe on the 
one hand and the rotating frame in fiat Minkowski spacetime on the other hand, differ 
substantially in their global behavior. 

This similarity can also be understood from the Taylor series expansion of (|4H) and (j44|) . 
Indeed, for small values of As or A'^ we find 



47rAi, 



V2 (^As - \aI + IIaI + 0(Ad) , (45) 
and 

47rA'|, = V2 (^A'^ - Ia% + U'^ + 0{A'^)^ . (46) 

Thus the two equations determining the Sagnac time delay only differ in the contribution 
of the fifth power of As or A'^. 
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Figure 4: Comparison between the Sagnac time delay for Godel's Universe (solid line) and 
for a rotating coordinate system in flat spacetime (dotted curve). Scaled proper times 
4:71 A\j and ^nA'^j of a light signal propagating from the origin to a mirror on the circle 
and back in Godel's Universe and in a rotating frame of reference in fiat spacetime, plotted 
versus the Sagnac time delays Ag and A'^ in scaled units. While the left ordinate applies 
to the absolute values of those scaled time delays, the dashed-dotted line magnifies the 
residual differences 4:n{A'^ — A|^) as indicated by the right ordinate. 

6 Summary and Outlook 

In conclusion we have investigated the Sagac time delay in Godel's Universe. Our analysis 
generalizes the results of our previous examination which were limited to small radial 
dimensions and rotation rates, to arbitrary sizes of the circular Sagnac interferometer and 
arbitrary rotation rates of Godel's Universe. The previous results exhibited a close relation 
between the Sagnac effect in Godel's Universe and in a rotating frame in fiat spacetime. 
By obtaining exact and invariant expressions for the Sagnac time delays in both systems, 
valid for arbitrary rotation rates and radii within G and G', we have demonstrated inhere, 
that this statement is also valid in general. 

The original Godel solution is not in agreement with cosmological observations of the 
background radiation jJSl- Therefore, one might be tempted to dismiss this remarkable 
exact solution of Einstein's field equations altogether. However, recently analogies be- 
tween light propagation in curved spacetimes and optics in moving media [21] or sonic 
propagation in condensed matter systems EH] have instigated renewed interest in this 
subject. For example, these articles consider experimental analogues of a gravitational 
black hole. Here, the event horizon emerges when the flow velocities of the optical medium 
or the condensate exceed the velocities of light or sound in these media, respectively. 
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How far can we push these analogies? Can we find an optical or sonic realization of the 
light propagation in Godel's Universe? Unfortunately, a global simulation of the Godel 
spacetime via an experimental analogue seems not to be very probable due to the causal 
structure of this metric discussed in Section I^Til But since we can approximate the local 
light propagation in Godel's Universe by a rotating coordinate frame in Minkowski space- 
time it should be possible to translate the physics of wave propagation in curved space 
locally to cold quantum gases. Vortices jSHl EZl HH] or alternatively Abrikosov lattices 
might offer a possible route to an analogue of Godel's Universe in the laboratory. 
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A Definitions and Conventions 

Throughout this article we use the signature (+, — , — , — ) of our metric g^^, with the 
determinant g = det g^i, and denote the covariant derivative with a semicolon and the 
ordinary partial derivative with a comma. 
The antisymmetric tensor 



V-9 

is defined in terms of the Levi-Cevita-Symbol 

1 for an even permutation 
^ ^ for an odd permutation ^ . (47) 

otherwise 

The components of the Riemann tensor 

— I p/^ pp pA* pP 

a/37 ~ "7./3 "/3.7 P/3 "7 P7 Cf/? 

result from the Christoffel symbols 

1 



2 

The Ricci tensor and the scalar curvature 



^ a/3 = o 3'^" i9ua,f3 + 5'i//3,Q " 5'a/3,i/) 



follow by contraction. 
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The covariant derivative of a contravariant vector field T" reads 



whereas for a covariant vector field it takes the form 

rp rp -pfj, rp 



B Integration of Null Geodesies 

In this Appendix we outline a procedure to obtain the solutions of the geodesic equations 
in Godel's metric. Since the main goal of the paper is the Sagnac time delay and in 
particular, an invariant formulation of it, we confine ourselves to special initial conditions 
appropriate for our experimental setup. 



B.l General Idea 

The central idea for solving the geodesic equations relies on finding simple expressions 
for the constants of motion. Some of these constants can easily be obtained if the metric 
possesses Killing vectors This feature can be understood by contracting the geodesic 
equations (|27|) with the contravariant components of a Killing vector, that is 

When we recall the Killing equation (j24|) and note, that the covariant derivative of a 
scalar is just the ordinary partial derivative, we find 



d 



Hence, every Killing vector corresponds to a constant of motion. 

In addition the norm of the tangent vector m^(A) is constant along the geodesies and 
therefore yields another constant of motion 

g^^u^u" = e\ (48) 

Since the metric coefficients fll8|l depend only on the radial coordinate r, the solution of 
the geodesic equations ()27j) can be found by making use of the first three Killing vectors 
in to obtain the simple constants of motion 

A = Mo(A) = c^u\\,) + r2(Ao) fiGw'(Ao) (49) 

B = u,{\) = r\Xo) Qgu^Xo) - r\\o) (^1 - (^^) j u\\o) (50) 

C = U3i\) = -u\Xo). (51) 
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Here Aq denotes the initial curve parameter and in the second step of these equations we 
have made use of the relation ■u^(Ao) = fi'^i/(Ao)w'^(Ao) with the metric coefficients (fTHj) . 
We substitute the expressions (jl^ - fjM^ for the covariant components of the four-velocity 



into the constant of motion 



and find with the contravariant metric coefficients 



/ 1 Hi-aT 



(9 










\ 



1 



(52) 



-1 / 



the equation 



^2 



(-)' 

\2aJ 



-(1 + 



+ 



G 



2AB 



\2a) 



ip \ 2, 

2a 



(53) 



Inserting Ui = QuU^ into this expression yields a differential equation in the radial coor- 
dinate r. The solution r(A) can then be used to formulate the differential equations for 
the coordinates t and by making use of the relations = g^^u^ and = g'^^u^. The 
solution of the third spatial coordinate z can be obtained from the last constant of motion 

m 

= g'^u-s = -C, z{X) = u%Xo){X - Ao) + ^(Ao). (54) 

Thus, the motion in z-direction is like a free motion in one dimension in flat Minkowski 
spacetime. 



B.2 Special Initial Condition 

In the remainder of this Appendix we focus on the solution of the geodesic equations for 
light, that is e = in (j48j) . subjected to the special initial conditions associated with the 
experimental arrangement considered in Section El In particular, the light pulse starts 
with the initial curve parameter Aq = at the position r(0) = 0, z{0) = and propagates 
in the z = plane. Therefore, we have to choose u^{0) = 0. We emphasize, that for r = 
the angle velocity m^(0) looses its meaning. Moreover, since we start at the origin and r 
is positive, we have to take u^(0) > 0. 

These special initial conditions reduce the constants of motion (pUjl - lf^ to 

A = uo{X) = cV(0) > (55) 
B = U2{X) = (56) 
C = ^i3(A) = 0. (57) 

These compact expressions allow us to integrate the differential equation for our 
radial coordinate r(A) as shown in the next section. 
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B.3 Radial Coordinate 

The constants of motion reduce the differential equation to 



which leads with ui = gn and u^{0) > to 



The two equations correspond to the two different signs of the radial velocity. 
Separation of variables 



r(A) 

dr 



±CM°(0) j dX 



|sin(ir/A)|. (59) 



V (2a) 
finally yields the radial geodesic 

r(A) 
2a 

Here we have introduced the periodicity rate 

rj = \/2u°(0)fiG (60) 

in such a way, that one cycle of the radial coordinate r(A) corresponds to the curve 
parameter A = 271/1]. 

It is apparent, that the solution is not differentiable at r = . This feature is due 
to the transformation to polar coordinates, which is singular at r = 0. Furthermore we 
recognize that these null geodesies touch the critical Godel radius tq = 2a for rjX = n. 



B.4 Coordinate Time 

The explicit expression (I59p for our radial coordinate allows us now to find the corre- 
sponding expressions for both the coordinate time t(A) and the angular coordinate 0(A). 
For this purpose we use the relation u° = g^^u^ with the contravariant metric coefficients 
(jH21) and the constants of motion and and arrive at 



„0(A) = n°(0) ^ \ '° i 
1 + ' 



o > 0. 

2a 



Hence, the coordinate time is increasing monotonously. 
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Substitution of (I^Uj) into the above equation leads to 



A 

t{\) = -u\0) A + 2u%0) [ (61) 

1 + sm [^r]X) 





which after integration yields 

2 
h 

Here we have used fjSU|) and introduced the integer 



t(A) = -M°(0) A + ^ (^arctan (^^2 tan (r7A/2) j + m(A) vrj . (62) 



m(A) 



r/A 1 
2^ 2 



(63) 



where the brackets denote the greatest integer less than or equal to the number inside 
them. The integer m guarantees the continuity, differentiability and monotony of t(A) for 
arbitrary A > 0. 

B.5 Angular Coordinate 

The dependence of the angular coordinate 0(A) on the curve parameter can be found 
in a way similar to the one of the coordinate time t{X). With the help of the relation 
^2 _ g'^f^ii^ and the constants of motion (|55|) and (|56p we arrive at 



u\X) 



Since the angular velocity v? looses its meaning at the point r = we perform the 
integration over one cycle only, that is < A < {2tt) /rj. In this case the integral 



0(A) = 0(0) +M°(0) U 



G 



dX 







1 + sin^ (ir/A) 

leads in analogy to under the restriction m(A) = 0, 1 to the angular coordinate 

0(A) = 0(0) + arctan (^^2 tan {r]X/2)^ + m(A) vr. 
We can rewrite this expression in order to allow arbitrary values of A > which results in 
0(A) = 0(0) + arctan (^^2 tan (r/A/2)) + {m{X) - m(^A - vr. (64) 
Equations (fS^. (jU^ and (jMj) represent the null geodesies for our special initial conditions. 
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